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Abstract 

We show that the Dual String Model with fusion leads, in heavy ion collisions, 
to strict saturation of the particle (pseudo-)rapidity density, normalised to the 
number of participant nucleons, as that number increases. Asymptotically, as 
y^ — > oo, with the number of participants fixed, this density approaches the 
nucleon-nucleon density. 

A comparison with recent WA98 data is presented. 



1 Introduction 

Recently, experimental [|l], H, 3 and theoretical |§, § papers have addressed the question 
of the dependence of measurable quantities (charged particle density, transverse energy, 
J/ip production rate) on the number A'^p^rt of participant nucleons in high energy heavy 
ion collisions. That information, as stressed in @|, as well as in ||5|, is extremely important 
as it allows for a better understanding of the initial conditions in the evolution of newly 
created dense matter and provides the information for discriminating among different 
models. 

In the analysis of the data of ||l| made in it was emphasised the role of the de- 
pendence of the particle rapidity density, dN/dy, normalised to the number of pairs of 
participants, and two models were considered. 

The HIJING Monte Carlo model ||^ contains soft and semi-hard interactions in an 
unitarised form, the pt threshold for jet production, po, being independent of energy and 
atomic mass number A. As in the model there is a growing number of hard collisions, 
growing faster than A^parti the quantity 



is an increasing function of A'^p^rt, or 



1 dN 
i^/'part dy 



d<l> 



^A/'part 



> 0. (2) 



The EKRT model 0, on the other hand, relates the hard physics threshold ps to a 
saturation criterion j^: the number of partons multiplied by their effective area (~ 1/Pt) 
must be less than the transverse area of interaction. The saturation transverse momentum 
Ps, in contrast to the parameter po of the previous model, is a growing function of ^/s and 
A. The assumption that physics is controlled by pt = Ps leads to the result that $ is a 
decreasing function of A^part , 

< 0. (3) 



^A/'part 



In Fig. |I|, taken from |||], we show the predictions for the two models, histogram for 
HIJING and dot-dashed line for EKRT. 

The Dual String Model (DSM) with string fusion was developed in |§. Concerning 
this model two remarks can be made: 

(i) There exist two components, a valence- valence component, proportional to A'^part, 
and a sea-sea (including gluons) component which grows faster than A^part and corresponds 
to multiple scattering. We have a structure somewhat similar to the HIJING model. 

1/3 

Without fusion DSM gives an increase of $ as, roughly, A'part- 

(a) With fusion, DSM leads to saturation of $ as A^part increases, not very different 
from what is obtained in the EKRT model (see, in particular, 0). 

In conclusion, in DSM with fusion, 

> 0, (4) 
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Figure 1: Central charged particle rapidity density per participating pair as a function of the number of 

participants. Results of HIJING (histograms), EKRT predictions (dot-dashed lines) and DSM predictions 

(solid lines) for central Au+Au collisions at ^/s = 56, 130, 200 AGeY. Also shown are results from pp and 

pp collisions and PHOBOS data (Everything in the figure except the DSM curves is taken from |Q|.) 



as in HIJING, and, as A'^part becomes larger and larger, 



^^part ^part- 



0, 



(5) 



The iVpart dependence of $ in DSM with fusion is also shown in Fig. |l| (solid lines). 
The parameters used are precisely those used in [^ to describe NA49 and PHOBOS data 

Ell. 



2 DSM and particle densities 



DSM is essentially the Dual Parton Model [^ with the inclusion of strings ||TT|. The 
strings may interact by fusing ||12| in the transverse plane of interaction and eventually 
one may reach a situation of percolation with the formation of extended regions of colour 
freedom, with the features of the expected Quark-Gluon Plasma. 

In the Dual Parton Model hadrons are considered as made up of constituents quarks 
(valence and sea quarks) and gluons. Two basic diagrams contribute to particle production 
(see 0). The valence- valence diagram corresponds to single inelastic scattering and the 
wounded nucleon model HTsf. The sea-sea diagram (including gluons) corresponds to 



the additional inelastic multiple scattering contributions. These contributions may be 
internal, parton multiple scattering within the original valence-valence contribution, or 
external, involving other nucleons. 

If h is the height of the valence-valence plateau, ah the height of the sea-sea plateau, 
2k the average number of strings produced in a nucleon-nucleon collision we can write 
1^, [l^ for the (pseudo-) rapidity particle density 



dN 
dy 



Na [2 + {2k - l)a] h + {p^ - NA)2kah, 



(6) 



NaNa 

where N^ is the number of pairs of participant nucleons 

1 



Na 



-N, 



part) 



(7) 



and ly^ is the average number of nucleon-nucleon collisions. From elementary multiple 



scattering arguments 14 we have 



Na 



N' 



4/3 



(8) 



One should notice that the number of nucleon-nucleon collisions is A^^ + {u^ 
and the number of strings is A''^ [2 + 2{k — 1)] + {u^ — NA)2k = 2ku^ . 
For Na = 1, nucleon-nucleon collision, Eq. @ gives 

= [2 + 2{k - l)a] h. 



-Na 



dN 
dy 



(9) 
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Assuming that h and a are energy independent (constant plateaus) the energy depen- 
dence of dN/dy\pp fixes the energy dependence of k. 
From Eq.s (|ID, (D, (0), (|) and (1) we can write 



HNa 



y) = 4>{V~s,y) + {N'f-l)2kah 



with 



0(V^,1/) 



dN 
dy 



(10) 
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The number of strings coming from nucleon multiple scattering — the second term in 
Eq. (H) — is Na{NJ — l)2k and they occupy the transverse interaction area S"^ , which, 
for central collisions, is approximately given by 



TT 1.14A^ 



1/3 



s. 



such that the dimensionless transverse density parameter r] is 



V 



1.14 



2kN^J^{NY^ 



1), 



(12) 



(13) 



where r^ ~ 0.2 fm is the string transverse section radius. Note that t] increases with A^i 
and \/s. 
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When fusion occurs, i.e., when strings cluster in the transverse plane, Eq. (|T0|) becomes 



^Na, V~s, y) = 0(v^, y) + F{7^) {Nf - l)2kah, 



where F{ri) is the particle production reduction factor |ff5 
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(14) 
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We shall next discuss the A^^, ^i and y dependence of $, Eq. (|T4|). 
For the Na dependence we have 
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and, making use of the definition of r], Eq. ( 
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As both (1 — e '') and {Nj — l)/Nj are increasing functions of Na, becoming constant 

as Na — > oo. 



9$ 



>0 



with 



5$ 

Wa 



as A^^ 



oo. 



(18) 



These are the results announced in the Introduction, Eq.s (H) and (||). 

The DSM with fusion thus predicts saturation of the particle rapidity densities per 
participant pair of nucleons, $, as Na increases. 

In Fig H we compare our model Eq. ([T^) with the very recent results (^ of WA98 
Collaboration on central Pb+Pb collisions at the CERN SPS. For large values of A'part 
the saturation is very clearly seen. Fig. ^ja, while for small values of A'part the rise of 
particle density with A'part, Fig. [^b, is well described. The values used for the parameters 
were h = 0.77 and a = 0.11, and for the pp density we used the same value as WA98, 



from |1^. Note that the parameters here have slightly different values in comparison to 
the ones used in 0. This is not surprising as the WA98 data are not compatible with 
NA49 data as given by PHOBOS H. 

Concerning the energy dependence of Eq. (|14D, if the energy is high enough 



^{Vs,y) = -f- 
dy 



2kah, 



(19) 
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and 



^NA,V~s,y)/(^{V~s,y) = l + F{r]){N'/ 



(20) 
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Figure 2: (a) Charged particle density per participant nucleon versus the number of participants; (b) 
absolute charged particle density versus the number of participants. The data from WA98 |2| refer to 
158 A GeV Pb+Pb collisions (filled circles), the open circle refers to pp collisions ||l9|; the solid line results 

from Eq. (|l|). 

As F{ri) goes to zero as ^/s increases we obtain the prediction that asymptotically the 
particle density per participant, for a fixed number of participants, approaches the nucleon- 
nucleon density and this approach is controlled by the function F{ri), Eq. ([T5|). This means 
that asymptotically heavy ion collisions, with respect to particle densities, become similar 
to nucleon-nucleon collisions. This is a consequence of the increasing role played by string 
fusion. 

In Fig. ^ we present our prediction for the dependence of the ratio (|20D on the energy 
for fixed values of A^part ■ 

Regarding the (pseudo-) rapidity dependence of $, Eq. ([14D, even without a more 
specific model, we can say that the valence-valence contribution dominates the large 
centre-of-mass \y\ region, as most of rapidity is taken by the valence diquarks. In the 
fragmentation region we expect the particle density per participant nucleon to be equal to 
the nucleon-nucleon density, 
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The sea-sea contributions are naturally shorter and more central in rapidity. 



3 Conclusions 

The DSM is a model with two components, the valence-valence component and the sea- 
sea component, the sea-sea component increasing its importance with energy and number 
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Figure 3: Ratio of the particle rapidity density per participating pair to the pp particle density, Eq. (| 
as function of cm. energy for fixed values of the number of participants, -/Vpart- 



of participants. This is somewhat similar to the HIJING Monte Carlo model, with soft 
and hard components. 

On the other hand, with fusion the DSM behaves, for large N a-, similarly to the EKRT 
model, but with strict saturation of the particle density per participant nucleon. However, 
in the original EKRT model the saturation criterion in the transverse plane is stronger 
than in case of fusion of strings. Here, saturation in the interaction area is asymptotic 
(when r] -^ oo) while in the EKRT model it occurs, using the string language, when 
r] = {nr^^) / {TTl.U'^N'^^){2kNAiN^^ - 1)) = 1. This causes the decrease of $ with Na in 
the EKRT original model. 

Probably different explanations, such as the ones based on string fusion, parton satu- 
ration, parton shadowing, are in some sense dual and refer to the same underlying physics 



14|| . What is becoming clear is that saturation of particle density puts strong constraints 



in models, and limits the rise of the (pseudo-)rapidity plateau at RHIC and LHC. 
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